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ABSTRACT 

We present the results from a study of the spatial distributions of cooling gas and in- 
trinsic X-ray absorbing material in a sample of nearby, X-ray bright cooling flow clus- 
ters observed with the Position Sensitive Proportional Counter (PSPC) on ROSAT. 
Our method of analysis employs X-ray colour profiles, formed from ratios of the surface 
brightness profiles of the clusters in selected energy bands, and an adapted version of 
the deprojection code of Fabian et al. (1981). We show that all of the cooling flow clus- 
ters in our sample exhibit significant central concentrations of cooling gas. At larger 
radii the clusters appear approximately isothermal. In detail, the spatial distributions 
and emissivity of the cooling material are shown to be in excellent agreement with 
the predictions from the deprojection code, and can be used to constrain the ages of 
the cooling fiows. The X-ray colour profiles also indicate substantial levels of intrinsic 
X-ray absorption in the clusters. The intrinsic absorption increases with decreasing 
radius, and is confined to the regions occupied by the cooling fiows. We explore a 
range of likely spatial distributions for the absorbing gas and discuss the complexities 
involved in the measurement of column densities from X-ray data. We show that the 
application of simple spectral models, in which the intrinsic absorber is treated as a 
uniform foreground screen, will naturally lead to significant underestimates of the true 
amounts of absorption. Comparison of our results with previously reported observa- 
tions made with Einstein Observatory Solid State Spectrometer (White et al. 1991) 
shows reasonable agreement, but requires that the absorbing material only partially 
covers the X-ray emitting regions. The masses of absorbing gas in the central (30 arc- 
sec radius) regions of the clusters, calculated under the assumption of solar metallicity 
in the absorbing material, can be accumulated by the cooling flows on timescales of 
a few 10* years, which are much less than the ages of the flows. This implies that 
most of the material deposited by the cooling flows in these regions cannot remain in 
X-ray absorbing gas. The results presented in this paper provide strong support for 
the standard model of inhomogeneous cooling flows in clusters of galaxies. 

Key v^rords: galaxies: clusters: general - cooling flows - intergalactic medium - 
X-rays: galaxies 



1 INTRODUCTION 

Clusters of galaxies are luminous X-ray sources. The bulk 
of the X-ray emission arises from bremsstrahlung and line 
radiation processes in the IntraCluster Medium (ICM). In 
the most luminous systems, the (2 — 10 keV) luminosity 
exceeds 10'*^ ergs~^ and the emission extends to radii > 2 
Mpc. The total mass of the ICM exceeds that of the galaxies 
by factors of 1 — 5. In the outer regions of clusters (r > 1 
Mpc) the ICM is highly diffuse (rie < 10"" cm"^) and the 
cooling time far exceeds the Hubble time (tcooi ^ lO^'^ yr). 
In the central regions of most clusters, however, the density 
of the ICM rises sharply and the cooling time inferred is 



significantly less than the Hubble time (fcooi ^ 10® yr; e.g. 
Edge, Stewart & Fabian 1992). In the absence of forces other 
than thermal pressure and gravity, the cooling of the ICM 
will lead to a slow net inflow of material towards the cluster 
centre; a process known as a cooling flow. (See Fabian 1994 
for a review of the theory and observations of cooling flows.) 

Within the idealized model of a homogeneous cooling 
flow (in which all of the gas at any particular radius has the 
same temperature and density), all of the cooling gas flows 
to the centre of the system where it is 'deposited', having 
radiated away its thermal energy. However, observations of 
clusters of galaxies show that cooling flows are not homo- 
geneous. The X-ray surface brightness profiles of clusters. 
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Table 1. Observation summary 



cluster 



OBSERVATIONS 
Date Exposure (s) 



X-RAY CENTROID (J2000) 
R.A. Dec. 



AbcU 85 


0.0518 


1992 


Jul 01 


10086 


00'M1™50, 


.5= 


-09° 17' 58' 


Abell 3112 


0.0746 


1992 


Dec 17 


7600 


03'' 17™ 58, 


.4= 


-44°14'27' 


Abell 426 


0.0183 


1992 


Feb 02 


4787 


Q3h]^gm48 


.5= 


41°30'27' 


Abell 478 


0.088 


1991 


Aug 31 


22139 


04'' 13^24, 


.9" 


10°28'04' 


Abell 496 


0.0320 


1991 


Mar 06 


8852 


04''33'"37, 


.9= 


-13°15'43' 


AbcU 644 


0.0704 


1993 


Apr 27 


10285 


08''17™25, 


.3= 


-06°49'59' 


Hydra A 


0.0522 


1992 


Nov 08 


18070 


09'n8™06, 


.2= 


-12°05'41' 


Coma 


0.0232 


1991 


Jun 17 


22108 


12''59™46, 


.1= 


27°56'21' 


Virgo 


0.0043 


1992 


Dec 17 


9961 


12''30™49, 


.8= 


12°23'32' 


Ccntaurus 


0.0104 


1994 


Jul 05 


3192 


12lM8™48, 


.7= 


-41°18'44' 


Abell 1795 


0.0634 


1992 


Jan 04 


36515 


13''48™53, 


.1= 


26°35'34' 


Abell 2029 


0.0767 


1992 


Aug 10 


12353 


15''10'"55, 


.7^ 


05°44'39' 


MKW3S 


0.043 


1992 


Aug 15 


9802 


15''21™51, 


.7" 


07°42'09' 


Abell 2199 


0.0309 


1993 


Jul 25 


40290 


16''28™38, 


.4= 


39°33'00' 


Cyg A 


0.057 


1993 


Oct 10 


9447 


19''59™28, 


.4= 


40°43'59' 


Scrsic 159 


0.0556 


1993 


May 28 


13137 


23''13™58, 


.3= 


-42°43'35' 


Abell 2597 


0.0824 


1991 


Nov 27 


7094 


23''25™20, 


.0= 


-12°07'31' 


Abell 4059 


0.0478 


1991 


Nov 21 


5514 


23''57™02, 


.0= 


-34°36'23' 



Notes: Summary of the PSPC observations. From left to right we list the names, rcdshifts, observation dates, efTective exposure times 
(after corrections for satellite dead time) and the coordinates of the peaks of the X-ray emission (in the 0.4 — 2.0 keV band) from the 
clusters. 



although sharply peaked, are not as peaked as they would 
be for homogeneous flows. Rather than all of the cool gas 
flowing to the centres of the clusters, the data show that 
cooled material is deposited throughout the central few tens 
to hundreds of kpc. Typically the deposition occurs with 
(r) oc r, where M{r) is the integrated mass deposition 
rate within radius r. Spatially-resolved X-ray spectroscopy 
of clusters also demonstrates the presence of a range of tem- 
perature and density phases in the central regions of cooling 
flow clusters {e.g. Fabian et al. 1994; Fukazawa et al. 1994; 
Allen, Fabian & Kneib 1996; Allen et al. 1996b; Fabian et al. 
1996). The X-ray spectra and imaging data firmly require 
that cooling flows axe inhomogeneous. 

For some years, the primary uncertainty with the stan- 
dard inhomogeneous model of cooling flows (Fabian 1994) 
was the fate of the cooled matter (see e.g. Fabian, Nulsen 
& Canizares 1991). No identification of the large masses of 
cooled gas, which should presumably be accumulated by the 
flows over their lifetimes, was made. However, a study of 21 
clusters observed with the Solid State Spectrometer (SSS) 
on the Einstein Observatory by White et al. (1991) showed 
that large masses of X-ray absorbing material are common 
in cooling flows. These authors also showed that signifi- 
cant absorption is not present in non cooling-flow clusters. 
Spatially-resolved X-ray spectroscopy of the massive cool- 
ing flow cluster Abell 478 with ROSAT (Allen et al. 1993) 
showed that the spatial distribution of the intrinsic absorb- 
ing material in this cluster is well-matched to that of the 
cooling flow. Allen et al. (1993) also demonstrated that the 
mass of absorbing material is in good agreement with the 
mass likely to have been accumulated by the cooling flow 
over its lifetime. The data thus suggest the X-ray absorbing 



material as a plausible sink for much of the cooled gas de- 
posited by cooling flows. Recent spectral studios of cooling 
flows with ASCA (Fabian et al. 1994; Fabian et al. 1996) 
have further confirmed the presence of cooling gas and in- 
trinsic absorbing material in cooling flows. However, to date, 
the constraints on the spatial distributions of these compo- 
nents remain poor. 

In this paper, we present detailed results on the dis- 
tributions of cool gas and intrinsic X-ray absorbing mate- 
rial in a sample of 17 cooling flow clusters observed with 
the Position Sensitive Proportional Counter (PSPC) on 
ROSAT. Our target clusters were selected on the basis of 
the following criteria; that they be X-ray bright (Fx, 2-10 > 
1.7x10"^^ erg cm~^ s~^ ; Edge et al. 1990), relatively nearby 
(2 < 0.1; allowing them to be resolved on scales of tens 
of kpc with the PSPC), and contain large cooling flows 
(M > 100 M© yr^^ ). We also include two lower- luminosity 
(smaller M) systems, the Ccntaurus cluster and the Virgo 
cluster, since they are previously well-studied and since good 
X-ray data exist for them. Finally, we included one non- 
cooling flow cluster, the Coma cluster, in our sample as a 
control. 

The large intrinsic column densities (AA'h ~ 10^^ atom 
cm~'^) in cooling flows inferred from the Einstein SSS (White 
et al. 1991) and ASCA data (Fabian et al. 1994) must have 
a significant effect on the flux detected with the PSPC at 
energies below about 0.8 keV. Indeed, few X-rays can es- 
cape the absorbed region at energies below the carbon edge 
in the detector window at 0.28 keV. However, this does not 
mean that no X-rays escape the cluster in the carbon band, 
since there is always a large contribution from the outer, 
less intrinsically absorbed, parts of the cluster along any 
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lino of sight. Given the spectral resolution of the PSPC 
[AE/E = 0.43(i;/0.93keV)"°-^] and the energy-dependent 
effects of absorption, one might expect that the lowest en- 
ergy range observable with the PSPC (0.1-0.4 keV) would 
be optimum for absorption studies. Unfortunately, however, 
the spatial resolution of the PSPC is substantially worse in 
this band than at higher energies. A significant fraction of 
the counts from a point source are scattered to radii > 1 ar- 
cmin in the 0.1-0.2 keV band, and in the 0.2-0.4 keV band 
the Point-Sprcad-Function (PSF) is still almost a factor of 
2 larger than at higher (0.4-2.0 keV) energies (Hasinger et 
al. 1992). In order to make meaningful spectra, the spatial 
regions selected must be at least as large as the largest rel- 
evant PSF. Since bright cooling flows typically have radii of 
only one or two arcmin, it is then difficult to make spatially- 
resolved studies of cooling flows with the PSPC using data 
from the softest band. Certainly, any study of the distri- 
bution of absorbing matter from such an exercise is very 
limited. 

The problems with the PSF are compounded by the 
large number of parameters required to model the spec- 
trum of an absorbed cooling flow, and the small number 
of emission-independent spectral bands of the PSPC. Even 
a simple model for an absorbed cooling flow involves 5 free 
parameters; the normalization, abundance and temperature 
of the ambient cluster gas, the cooling rate in the cluster 
and the column density of intrinsic absorption. (This as- 
sumes that cluster redshift and Galactic column density are 
known beforehand.) Although the counts in each of the 256 
PSPC spectral channels are statistically independent, the 
limited spectral resolution of the detector means that there 
are only about 5 independent spectral bands from a mod- 
elling point of view. Consequently, most spectral analyses 
with ROSAT data employ simple models of gas at a single 
temperature and with a uniform screen of absorption. We 
term these 'single- phase' models. [Temperatures and column 
densities determined with single-phase models should only 
be interpreted as some emission/detector- weighted mean 
over a range of temperatures or column densities at any 
projected radius.] Only with the much higher spectral reso- 
lution of instruments such as the Einstein Observatory SSS 
and the ASCA Solid-state Imaging Spectrometers (SIS) can 
more detailed spectral fits be attempted. 

In this paper we make use of X-ray colour profiles, 
rather than individual spectra, to study the spatial distri- 
butions of cool gas and absorbing material in clusters. By 
restricting our study to the 0.4-2.0 keV band, where the 
FWHM of the PSF is approximately constant at 25 arc- 
sec, we can make optimum use of the spatial resolution 
of the PSPC. Emission lines, principally due to Fe-L tran- 
sitions characteristic of gas cooling through temperatures 
of ~ 1 keV, are strongest between energies of 0.8 and 
1.4 keV. Intrinsic column densities at the level of ~ 10^^ 
atom cm~^ should have a marked effect on the spectrum 
below about 0.8 keV. The remaining 1.4-2.0 keV band is 
comparatively little-affected by emission lines and absorp- 
tion. We have therefore divided the 0.4-2.0 keV band into 
bands B, covering the range 0.41-0.79 keV, C from 0.80- 
1.39 kcV and D from 1.4-2.0 kcV. The ratio of B to D should 
be sensitive to absorption, and the ratio of C to D to the 
presence of cool gas (where 'cool' implies a temperature sig- 
nificantly less than that of the bulk of the ICM). By creat- 
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Figure 1. Representation of the energy bands used in the X-ray 
colour profile analysis, overlaid on the effective area curve for the 
PSPC. 

ing ratios of the surface brightness profiles of the clusters in 
these bands, we parameterize the amounts of cool gas and 
absorbing material in the clusters and, to the limits of the 
spatial resolution of the PSPC, constrain the spatial distri- 
butions of these components. We then compare these results 
to the predictions from standard cooling flow models. 

The structure of the paper is as follows. In Section 2 we 
summarize the observations and data reduction. In Section 3 
wo present the results from a standard doprojection analysis 
of the data. In Section 4 we present the results on the X-ray 
colour profiles and relate them to expectations from simu- 
lated cluster spectra. In Section 5 we compare the observed 
colour profiles to the predictions from the doprojection anal- 
ysis and determine the results on the distributions of cool 
gas and X-ray absorbing material in the clusters. In Section 
6 we discuss the complexities associated with the measure- 
ment of intrinsic column densities with PSPC data. Section 7 
summarizes our conclusions. Throughout this paper, _ffo=50 
kms~^ Mpc~^, Q = 1 and A = are assumed. 



2 OBSERVATIONS 

The PSPC observations were carried out between 1991 
March 6 and 1994 July 5. Exposure times range from 3192s 
for the Centaurus cluster to 40290s for AbeU 2199. Details of 
the individual cluster observations are summarized in Table 
1. In cases where more than one observation of a source was 
made, only the longest of the individual exposures has been 
used. This minimizes systematic errors in the central surface 
brightness profiles that could be introduced by mosaicing. 
For the Coma cluster, for which four individual observations 
were made at different pointing positions, we have used the 
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Table 2. Definition of the energy bands 



Band Energy range (keV) 



Ai 0.10 - 0.40 

A2 0.20 - 0.40 

B 0.41 - 0.79 

C 0.80 - 1.39 

D 1.40 - 2.00 

F 0.41 - 2.00 



data set that places the peak of the cluster emission closest 
to the centre of the PSPC field of view. For the Centaurus 
cluster, we have used only the data from the 3.2ks observa- 
tion in 1994 July, which was less-affoctcd by scattered solar 
X-ray emission than earlier observations. 

The data were reduced using the Starlink ASTERIX 
package. X-ray images with 15 x 15 arcsec^ pixels were ex- 
tracted in each of the energy bands listed in Table 2. Regions 
affected by emission from point sources were masked out and 
ignored. Corrections for telescope vignetting were made and 
the data were background subtracted using circular regions 
of radius 6 arcmin, offset by ~ 40 arcmin from the centres 
of the fields (except in the case of the Virgo cluster, where a 
3 arcmin radius region offset by 45 arcmin was used) . Accu- 
rate centroids for the X-ray emission were determined from 
the full 0.4-2.0 keV (Band F) images. These are also listed 
in Table 1. X-ray surface brightness profiles were then ex- 
tracted, about these centroids, in each energy band. 



3 DEPROJECTION ANALYSIS 

The F band (0.4-2.0 keV) images of the clusters have been 
analysed using an updated version of the doprojcction code 
of Fabian et al. (1981; see also Thomas, Fabian & Nulsen 
1987). Under assumptions of spherical symmetry and hydro- 
static equilibrium in the ICM, the deprojection technique 
can be used to study the properties of the intracluster gas 
{e.g. density, pressure, temperature, cooling rate) as a func- 
tion of radius. The deprojection method requires that ei- 
ther the total mass profile (which defines the pressure pro- 
file) or the gas temperature profile be specified. Following 
ASCA observations of a number of the clusters included in 
the sample (Fabian et al. 1996; Ohashi et al. 1996), and 
the results from the combined X-ray and gravitational lens- 
ing study of the cooling-flow cluster PKS0745-191 by Allen 
et al. (1996a), we assume that the mass- weighted tempera- 
ture profiles in the clusters remain approximately constant, 
at the temperatures listed in Table 3. Note that the as- 
sumption of a constant mass- weighted temperature profile is 
consistent with ROSAT results which often show a drop in 
the emission-weighted temperatures in the central regions of 
nearby cooling-flow clusters. Such measurements will natu- 
rally result from the presence of high emissivity cooling gas, 
distributed throughout the (inhomogeneous) cooling flow 
with M 25 r, even when the mass- weighted temperature pro- 



Table 4. A comparison of isothermal and cooling gas 



Isothermal C. Flow 

kT (keV) C/D B/D C/D B/D 



1 


5.30 


2.08 


7.76 


4.76 


2 


3.27 


1.61 


4.72 


2.60 


3 


2.97 


1.51 


4.02 


2.21 


4 


2.83 


1.44 


3.71 


2.02 


5 


2.74 


1.39 


3.53 


1.88 


6 


2.68 


1.36 


3.40 


1.79 


7 


2.64 


1.34 


3.30 


1.73 


8 


2.60 


1.32 


3.22 


1.68 


9 


2.58 


1.30 


3.15 


1.66 


10 


2.56 


1.29 


3.11 


1.63 



Notes: A comparison of the theoretical C/D and B/D ratios 
for isothermal gas and a constant pressure cooling flow. For 
the cooling flow data, the kT value is the ambient cluster tem- 
perature from which the gas cools. A Galactic column density 
of 10-^'' atom cm~-^ and a metallicity of 0.5Zq are assumed. 

file of the infiowing material remains constant {e.g. Allen et 
al. 1996a,b). 

The primary results on the cooling flows in the clus- 
ters; the cooling times in the central 30 arcscc bin, the cool- 
ing radii (the radii at which the cooling time first exceeds 
the Hubble time) and the integrated mass deposition rates 
within the cooling radii, are listed in Table 3. In all cases 
the mass deposition from the cooling flows is distributed 
throughout the cooling radii of the clusters, with M 2$ r. 



4 THE X-RAY COLOUR PROFILES 

As discussed in Section 1, our method of analysis makes 
use of X-ray colour profiles, rather than the more standard 
approach of taking annular or regional spectra, to study the 
distributions of cool gas and intrinsic absorbing material in 
the clusters. Such a technique allows us to make the best 
use of the spatial resolution offered by the ROSAT PSPC. 

We have constructed radial surface brightness profiles 
for the clusters in each the energy bands, B, C and D (Table 
2). The surface brightness profiles were binned to a spatial 
scale of 30 arcsec pixel" (slightly larger than the FWHM 
resolution of the PSPC) which provides a good signal-to- 
noisc ratio in the pixels over the regions of interest. Figs. 2 
and 3 summarize the results on the X-ray colour profiles. Fig. 
2 shows the C/D profiles, which are approximately constant 
at large radii but increase in their central regions. The B/D 
profiles, shown in Fig. 3, are generally flatter than the C/D 
data. 

Figs. 4 and 5 illustrate how the X-ray colour profiles 
may be used to probe the physical properties of the clusters. 
These figures show the behaviour of the C/D and B/D ratios 
in simulated cluster spectra as a function of temperature, 
column density and metallicity. Results arc shown both for 
isothermal gas and for gas cooling at constant pressure, from 
the ambient cluster temperature to a temperature below the 
ROSAT band i. e. a constant pressure cooling fiow. The sim- 
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Table 3. Cooling flow results 





kT 




''cool 


^ cool 


^cool 


M 




(keV) 


(10^0 atom cm-2) 


(kpc) 


(arcmin) 


(109 yr) 


(Moyr-l 


Abell 85 




2.96 


137+2 , 

-34 


1 66+2- 
""-0.41 


2 85+°-;? 

""-0.11 


187+',? 

-42 


Abell 3112 


4 1+2-3 


2.78 


251+^1 


^•^"-0.40 


2 l3+°-°7 

^•^'^-0.08 


406+'^^ 


Abell 426 


"■■^-0.2 


13.7 


169+^ 


5 46+°-?? 


1 06+°-°? 


456+;? 

— 19 


Abell 478 


.^0+0.7 


20.0 


213+1'' 


1 61+"-!'' 

-^■"-^-0.36 


2 59+°-™' 


' ^ ' -107 


Abell 496 


4 8+°-t 

^•°-Q.8 


4.27 


105+" 


1 99+°-?6 

^•""—0.26 


1 84+?-°2 

^•"^—0.09 

8 28+i-?„ 

"■''"— 0.80 


85+i° 


Abell 644 


8.1+2-6 


8.61 


200+^? 

•'""— 61 


1 83+5-t5 


244+63 


Hydra A 


3 6+°'? 

''■"-0.5 


4.84 


161+?? 


1 93+S-?i 


1 63+°-'?! 

i.uo_„ Q4 


267+i*? 


Coma 


8.0+0-3 

—0.3 


0.90 


otf 


0+0.7B 


18.4+?\,i 

— 5.9 


— 


Virgo 


2 4+0.3 


1.71 


98+2* 


13.16lf-? 


21+°-°i 


37lg** 


Centaurus 




8.87 




4.88l?-« 


70+°-°6 

'^•"^-0.04 


27+^ 


Abell 1795 


c 1+0.4 
"■^-0.5 


1.11 


1921?^ 




2 14+0.05 


447+50 


Abell 2029 


7 1+2.0 
'■^-1.4 


3.05 


2001^3 


1 7o+°-05 


^•'"-0.08 


SSOlg^ 


MKW3s 


9 «+0-7 


2.86 


15411 


r, r,i+1.05 

^•^1-0.46 


9 17+0.15 
^•-^'-0.17 


161+^^ 


Abell 2199 


4 7+0-4 
^•'-0.4 


0.87 


1421^4 


9 70+0.47 
^•'^-0.53 


1 94+°-°4 
^•^^-0.04 


162+26 


Cygnus A 


^■'--1.3 


32.8 


157116 


1 74+0.51 
^- '^-0.49 


1 39+°-°5 


32OIS 


Sersic 159 


2 9+11 


1.80 


167lf? 


1 89+°-^'^ 
-^■''^-0.12 


1 7C+0.06 


23ii;i 


Abell 2597 


> 3.8 


2.46 


162+^'^ 


1 29+°-*6 


2 44+°-" 

^•^^-0.07 


299+?7 


Abell 4059 


q c+0.6 


1.61 


134+" 


1 74+1-00 
^- '^-0.49 


q crq+0.56 


1151^^ 



Notes: A summary of the results from the dcprojection analyses of the F band images. Temperature (kT) data arc from Edge et al (1990) 
except for Abell 426 (Allen et al. 1992) and Abell 478 (Allen et al, 1993). Nn values are Galactic column densities from Stark et al. 
(1992) except for Abell 478 which is from Johnstone et al. (1992). Cooling radii (r-cool) ^fe quoted in both kpc and arcmin. Cooling times 
(tcool) for the central (30 arcsec radius) bin are in 10^ yr, and mass deposition rates (M ) in Mq yr~^ . Errors on the cooling times are 
the 10 and 90 percentile values from 100 Monte Cairlo simulations. The upper and lower confidence limits on the cooling radii are the 
points where the 10 and 90 percentiles exceed and become less than the Hubble time, respectively. Errors on the mass deposition rates 
are the 90 and 10 percentile -values at the upper and lower limits for the cooling radius. 



ulated spectra were constructed using the plasma code of 
Kaastra & Mewe (1993) and the photoelectric absorption 
models of Morrison & McCammon (1983). Wc adopt the def- 
inition of solar abundance from Anders & Grevesse (1989). 
In all cases the absorbing column density is modelled as a 
uniform screen, with solar metallicity, at zero redshift. A 
cluster redshift of 2 = 0.01 is assumed (although the results 
do not vary significantly for any redshift < 0.1). 

The theoretical curves for an isothermal plasma are 
shown in Figs. 4(a) (d). At high temperatures {kT > 3 keV) 
both the C/D and B/D ratios arc weak functions of temper- 
ature. At lower temperatures, however, the C/D ratio (in 
particular) increases sharply as the temperature drops, due 
primarily to Fe-L emission lines in the C band produced by 
gas at ~ 1 keV. Both the B/D and C/D ratios decrease with 
increasing column density, although the variation of the B/D 
ratio is steeper (particularly at low column densities). Both 
ratios also exhibit some variation with metallicity. In gen- 
eral, as the metallicity rises the C/D ratio increases (slightly) 
and the B/D ratio drops, although for kT > 3 keV the func- 
tions are very weak. 

The behaviour for cooling gas, shown in Figs. 5(a)-(d), 
is qualitatively much like that for isothermal emission, ex- 
cept that the C/D and B/D ratios generally give higher val- 
ues for the same temperature, column density and metal- 
licity. This is illustrated by the data presented in Table 4 
where we list the C/D and B/D values for isothermal and 



cooling plasmas as a function of temperature, for a fixed col- 
umn density of 10^° atom cm"^ and a metallicity of O.5Z0. 
The variation of B/D with metallicity, Fig. 5(d), is fiatter 
for cooling gas than for an isothermal plasma, and exhibits a 
downwards trend at low metallicities not seen in the isother- 
mal data. 

In what follows, we use the C/D profiles to demonstrate 
the presence of distributed cooling gas in the central regions 
of the cooling fiow clusters, and the B/D ratio to demon- 
strate the presence of intrinsic X-ray absorbing gas. 



5 RESULTS FROM THE COLOUR PROFILES 
5.1 C/D and cooling radii 

The C/D colour profiles shown in Fig. 2 exhibit a num- 
ber of striking features. Firstly, we see that at large radii 

{1' > '"cool) the profiles arc essentially fiat. The observed 
values arc generally in good agreement with the theoretical 
predictions for isothermal gas at the temperatures listed in 
Table 3, absorbed by the appropriate Galactic column den- 
sity (although significant deviations are observed for Abell 
478, the Coma cluster, and Abell 2199, the origins of which 
are discussed in Section 5.3.) The observed and predicted 
C/D and B/D ratios at large radii are summarized in Table 
5. 



© 0000 RAS, MNRAS 000, 000-000 



6 S.W. Allen and A.C. Fabian 





Abell 496 



Abell 644 




2 4 
Radius (arcmin) 



Q 
\ 
O 




1 2 3 
Radius (arcmin) 



Figure 2. The C/D (0.80 — 1.39 kcV/1.40 — 2.00kcV) X-ray colour profiles (with Icr error bars). The cooling radii determined from the 
dcprojection analysis (Section 3) are indicated by the vertical dot-dashed lines. The 10 and 90 per cent confidence ranges on those radii 
are marked by the dotted lines. The profiles are essentially flat at large radii but show clear increases within r^ool- 



Within radii r < rcooi, however, the cooling-flow clus- 
ters exhibit clear increases in their C/D ratios. From com- 
parison with the theoretical curves in Figs. 4 and 5, we see 
that these increases can only plausibly be explained by the 
presence of cooler gas. The X-ray colour profiles thus pro- 
vide firm evidence for distributed cool(ing) gas in the central 
regions of the cooling flow clusters and show that this gas 
is confined within the cooling radii inferred from the de- 
projection analyses. For the one non-cooling flow cluster in 



the sample, the Coma cluster, no significant enhancement 
in the C/D ratio (outside the central 30 arcsec i.e. 20 kpc) 
is observed. 



5.2 The spatial distribution of cooling gas 

The C/D profiles shown in Fig. 2 clearly demonstrate the 
presence of distributed cool gas in the centres of the cool- 
ing flow clusters. The regions over which the cool gas is 
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Figure 2. - continued 



distributed are well matched to the cooling radii inferred 

from the deprojcction analysis (Section 3). Given this agree- 
ment, we have attempted a more sophisticated comparison 
of the C/D profiles with the predictions from the deprojcc- 
tion code. 



The deprojcction analysis employed in Section 3 divides 
the X-ray emission from the cluster into a series of concentric 
spherical shells. The luminosity in shell j may be written as 
the sum of four components. 



Li = AMiHi + AM, A$,- -|- 



(1) 



where AMj is the mass deposited in shell j, Hj is the en- 
thalpy of the gas in shell j, and A$j is the gravitational 
energy released in crossing the shell. X^^Zj AM, is the mass 
flow rate through shell j, and AHj the change of enthalpy 
of the gas as it moves through the shell. 

Thus the first term in equation 1 accounts for the en- 
thalpy of the gas deposited in shell j. The second term is the 
gravitational work done on the gas deposited in shell j. The 
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Figure 2. - continued 



third and fourth terms account for the gravitational work 
done on gas flowing through shell j to interior radii, and the 
enthalpy released by that gas as it passes through the shell. 
Since we have constrained the mass-weighted temperatures 
in the cooling flows to remain approximately constant with 
radius (at the values listed in Table 3) the luminosity in the 
fourth term is negligible. 

The luminosity contributed by the second and third 
terms will have a spectrum that can be approximated by 
an isothermal plasma at the appropriate temperature for 
the cluster (c/. Fig. 4). The emission accounted for in the 
first term, however, should have the spectrum appropriate 



for gas cooling from the ambient cluster temperature (c/. 
Fig. 5). By assigning the appropriate C/D values from the 
theoretical curves, corrected for redshift effects, to each of 
the terms in equation 1, the deprojection code can be used 
to predict intrinsic and projected C/D colour profiles for 
the clusters. These predictions can then be compared to the 
observed values. 

In Fig. 6 we show the comparison of the observed and 
predicted C/D profiles for the cooling flows in the sample. 
Due to the finite ages of the clusters, we limit the compari- 
son to those radii where the cooling time of the cluster gas 
is < 5 X l(f yr. In the outer regions of the clusters the emis- 
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Figure 3. The B/D (0.40 — 0.79 kcV/0.80 — 1.39 kcV) X-ray colour profiles (with Icr error bars). Circled points mark the values predicted 
from the deprojcction code (Section 5.3) assuming Galactic absorirtion. These jHcdictcd values arc only plotted for radii where the cooling 
time of the gas is less than 5 X 10^ yr. Note the clear deficit in the observed B/D values relative to the predictions from the deprojection 
code in the central regions, which indicates the need for excess absorption. In the outer regions the circled points mark the best-fitting 
constant B/D value (as listed in Table 5). 



sion is assumed to be isothermal, at the temperature given 
in Table 3, and absorbed by the appropriate Galactic col- 
umn density. The predicted curves are then re-normalized 
to match the observed C/D values in the outer regions of 
the clusters. (Note that the normalization factors, which arc 
given by the ratios of the observed to predicted C/D ra- 
tios in the outer regions of the clusters, are typically close 
to unity.) We do not attempt any comparison for either the 



Centaurus or Virgo clusters since both systems are relatively 
cool and are known to contain strong metallicity gradients 
(Allen & Fabian 1994; Fukazawa et al. 1994). The low red- 
shifts of these systems also means that some drop in their 
mass-weighted temperatures, as the gravitational potentials 
of the clusters merge with those of the central galaxies, are 
likely in the innermost 1-2 bins. Indeed, such a temperature 
drop is required to explain the large C/D ratio observed in 
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Figure 3. - continued 



the central bin of the Centaurus data. We therefore defer de- 
tailed study of the data for the Virgo and Centaurus clusters 
to future work. (The relatively poor quality of the data for 
Abell 3112 and Abell 4059 also prohibit any detailed com- 
parison with the deprojection code.) We note that hotter 
clusters do not generally exhibit metallicity gradients {e.g. 
Ohashi et al. 1996) and the effects of metallicity gradients 
on the X-ray colour profiles of such systems will, in any case, 
be small (cf. Figs. 4,5). 

The agreement between the observed and predicted 
C/D ratios is remarkably good. For the best data sets (Abell 



426, Abell 1795 and Abell 2199; see also Fig. 13) the agree- 
ment is especially good. In all other cases the predicted C/D 
values either agree with, or slightly under-predict, the ob- 
served values. Under-prediction in the central bin may again 
indicate a small drop in the ambient cluster temperature. 

5.3 The B/D ratios and intrinsic absorption 

As with the C/D results, the B/D colour profiles shown in 
Fig. 3 are essentially flat at large radii (r > rcooi)- The ob- 
served and predicted values at large radii are generally in 
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Table 5. C/D and B/D values at large radii 









C/D 




B/D 




cluster 


radii 


observed 


predicted 


observed 


predict* 




3.0- 


-6.0 


9 c;« 4- n 

^.OO HZ u.uo 


2.58 


X.^O ZIZ u.uu 


1.20 


.nut:;!! OXJ.^ 


3.0 - 


-6.0 




2.80 


1.34 ± 0.11 


1.32 


Ahpll 42fi 


6.0 - 


- 15.0 


9 ■11 + n 03 


2.24 


n «2 + n 01 


0.74 


Ahpll 478 


2.5 - 


-4.5 


-L.oo ^ yj.\jo 


2.07 


n QQ _|_ n 09 


0.58 




3.6 - 


- 9.0 




2 64 


1 1 Q _!_ n r\A 

L.IO ^ U-U4: 


1 1 7 

J — L / 


Aholl (\Ad 


3.0 - 


- 7.2 




2.35 


u.oy zjz u.ui 


0.90 


xiycira a 


3.0 - 


- 4.5 


z. i z it U. iU 


Z. t 1 


l.ZD it U.Uo 


1 on 


Coma 


2.0 - 


12.0 


2.69 ±0.02 


2.61 


1.42 ±0.01 


1.33 


Virgo 


10.0 - 


- 12.0 


3.04 ± 0.05 


3.07 


1.57 ±0.04 


1.48 


Centaurus 


6.0- 


■ 12.0 


2.66 ± 0.09 


2.58 


1.07 ± 0.04 


0.97 


Abell 1795 


2.5 - 


-5.5 


2.67 ±0.04 


2.76 


1.46 ±0.03 


1.40 


Abell 2029 


3.0- 


-4.5 


2.46 ± 0.09 


2.58 


1.19 ±0.05 


1.19 


MKW3S 


3.0- 
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2.85 ±0.19 
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Abell 2199 
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1.74 ±0.11 


1.48 
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1.8 - 


- 3.0 


2.79 ±0.21 


2.83 


1.34 ±0.16 


1.35 
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4.2 - 


-8.4 


2.79 ±0.17 


2.91 


1.26 ±0.09 


1.44 



Notes: The observed and predicted C/D and B/D ratios in the outer regions of the clusters. Column 2 lists the range of radii (in arcmin) 
that define the 'outer regions'. Columns 3 and 5 summarize the observed C/D and B/D values and la errors. Columns 4 and 6 list 
the predicted ratios for an isothermal plasma, at the temperature and redshift of the cluster. Galactic column densities (Table 3) and 
a metallicity of O.5Z0 are assumed. The general agreement between the observed and predicted C/D and B/D ratios demonstrates the 
validity of the isothermal assumption in these regions, and the absence of excess absorbing material at these radii (see discussion in 
Section 5.3). 



Table 6. Absorption results from the B/D ratio 



Cluster 






AArH(i) 


AiVH(2) 


Abell 85 


2.96 


2.34 ± 0.60 


5 4+1-'^ 


10.5 ±2.7 


Abell 426* 


13.7 


0.50 ± 0.55 




2.3 ± 2.5 


Abell 478 


20.0 


7.4 ±0.63 


20.1111 


33.3 ± 2.8 


Abell 496 


4.27 


2.57 ±0.62 




11.6 ± 2.8 


Hydra A 
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3.07 ±0.41 


7 5+1-"' 
'•^-1.2 


13.8 ± 1.8 


Abell 1795 


1.11 


2.20 ± 0.20 


5 0+°-^ 
^■^-0.5 


9.9 ±0.9 


Abell 2029 


3.05 


0.83 ±0.42 


1 7-1-1.0 

-0.9 


3.7± 1.9 


MKW3s 


2.86 


7.29 ±0.88 


11 7+1-8 


32.8 ±4.0 


Abell 2199 


0.87 


2.67 ±0.25 


6 3+°-* 
"■•^-0.7 


12.0 ± 1.1 


Cygnus A 


32.8 


5.5 ± 1.0 


1 ■? 4+3-3 

10.1_2g 


24.8 ± 4.5 
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1.80 


1.87 ±0.48 


4 2+1-2 


8.4 ±2.2 


Abell 2597 


2.46 


2.36 ± 0.58 




10.6 ±2.6 



The intrinsic column densities and Icr errors (in units of 10-^" atom cm^-^) for the central 30 arcsec regions of the clusters, inferred from 
the B/D analysis (Section 5.3). Column 2 lists the Galactic column densities. Column 3 (AN^) summarizes the excess column densities 
determined with the simple uniform absorbing screen model. Column 4 [AAfH(l)] lists the corrected intrinsic column densities, determined 
from equation 4, with the partial covering model (the preferred model). A covering fraction, / = 0.5, and a photon energy, E = 0.4 keV, 
have been assumed for all clusters except Abell 478 and Cygnus A (which have much higher Galactic column densities) for which E = 0.5 
keV is adopted. For MKW3s the observed AAT^ is inconsistent with / = 0.5 and / = 0.75 is used instead. Column 5 [AArH(2)] lists the 
corrected intrinsic column densities determined with equation 5, for the multilayer absorption model. An intermediate correction factor 
of 4.5 has been assumed. * Between radii of 1-2 arcmin in Abell 426, the corrected intrinsic column density is AiVH(l) ~ 5 — 10 x 10^" 
atom cm~2. 
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Figure 3. - continued 



good agreement (although with a few exceptions detailed 

below) indicating that the omission from these regions can 
be well-described by a simple isothermal model, with the 
temperatures and Galactic column densities listed in Table 
3. At small radii, many of the B/D profiles increase. How- 
ever, the gradients are not as steep as for the C/D results, 
and for a number of the clusters the B/D profiles are es- 
sentially flat or even decrease slightly in the innermost bins 
{e.g. Abell 478). 

The most notable exception to the agreement between 
the observed and predicted colour ratios at large radii is 



Abell 478, for which the B/D (and C/D) ratio implies a 
Galactic column density of --^ 3 x 10'^^ atom cm~'^, an excess 
of ~ 10^^ atom cm ^ over the value listed in Table 3. In 
contrast, the data for Abell 426 require a column density 
of ~ 1.1 X 10^^ atom cm , slightly less than the Galactic 
value. For the Coma cluster, the B/D and C/D ratios imply 
an emission-weighted temperature at largo radii of 6 kcV, 
somewhat less than the 8 keV listed in Table 3. For Abell 
2199 the data suggest an emission weighted temperature of 
~ 3.5 keV and a column density of ~ 5 x 10^^ atom cm~^. 
Similarly, the data for the Virgo cluster imply a temperature 
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Figure 4. Theoretical predictions for the C/D and B/D ratios for isothermal gas at temperatures of 1-10 keV. (a) C/D as a function 
of column density {Z = 0.5Zq fixed), (b) B/D as a function of column density {Z = 0.5Zq fixed), (c) C/D as a function of metallicity 
(ATjj = 10^" atom cm~^ fixed), (d) B/D as a function of metallicity (ATjj = 10^" atom cm~^ fixed). In all cases the temperature of the 
gas increases with decreasing C/D and B/D along the y axis. 



of ~ 2 keV and a column density of ~ 1.3 x 10^° atom cm~^, 
both slightly less than the nominal values. 

Following our success in applying the dcprojection code 
to the prediction of the C/D profiles (Section 5.2), we have 
used the same technique to predict the B/D profiles. The 
predicted data are shown as the circled points in Fig. 3. The 
results demonstrate a clear tendency for the predicted values 
to overestimate the observations in the central regions. 

The theoretical curves presented in Figs. 4 and 5 show 
that the most plausible interpretation for the central deficits 
in the B/D observations, with respect to the predicted data, 
is the presence of intrinsic absorbing material in the clusters. 
The excess column densities required to bring the observed 
and predicted B/D values for the central bins into agreement 
are summarized in Table 6. These values have been calcu- 



lated using XSPEC models for the projected spectrum of 
the central 30 arcsec (radius) in each cluster, as implied by 
the deprojection analysis. The excess absorption has then 
been modelled as a uniform screen in front of the cluster. 
Note that since the observed C/D ratios for the central bins 
sometimes (slightly) exceed the predicted values, the column 
densities listed in Table 6 may (slightly) underestimate the 
true values. 

Where comparisons between the observed and predicted 
B/D ratios are possible for more than just the central bin, 
the data suggest that the column density increases with de- 
creasing radius within the cooling flow. This is consistent 
with a picture wherein the intrinsic absorber is accumulated 
from material deposited by the cooling flow. The only ex- 
ception is Abell 426 (the Perseus Cluster) for which the ap- 
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Figure 5. Theoretical predictions for the C/D and B/D ratios for gas coohng at constant pressure from and upper temperature of 1-10 
keV to a lower temperature of .001 keV. (a) C/D as a function of column density {Z = O.BZq fixed), (b) B/D as a function of column 
density {Z = 0.5Z© fixed), (c) C/D as a function of metallicity {N^i = 10^" atom cm~^ fixed), (d) B/D as a function of metallicity 
(ATjj = 10^" atom cm~^ fixed). In all cases the temperature of the gas increases with decreasing C/D and B/D along the y axis. 



parent excess column density at radii of 1-2 axcmin is 5 — 10 
times larger than for the central bin. It is likely that the re- 
sults for the central regions of Abell 426 are affected by the 
complex morphology of the cluster (Bohringer et al. 1993) 
and emission from the active nucleus of the central galaxy. 
Note also that although the presence of temperature and 
metallicity gradients complicate the interpretation of the 
B/D data for the Virgo and Centaurus clusters, the com- 
parison of their B/D ratios, which are only observed to rise 
within the central arcminute or so, with their C/D ratios, 
which in contrast rise steadily within radii of 7 and 3 arcmin 
respectively, again indicates the presence of distributed ab- 
sorbing material. 

Our results confirm that excess absorption is common 
in cooling flows. All of the cooling flow clusters in our sam- 



ple, for which an X-ray colour profile/deprojection compar- 
ison is possible, show evidence for intrinsic absorption. Note 
also that the cooling flow with the lowest intrinsic column 
density, Abell 2029, is also the only cluster with a Central 
Cluster Galaxy (CCG) that does not exhibit optical emis- 
sion lines [e.g. Hu, Cowie & Wang 1985; Johnstone, Fabian 
& Nulsen 1987). This suggests a link between the presence 
of X-ray absorbing gas in cooling flows and the production 
of optical emission lines in CCGs. Such a result is compati- 
ble with models in which the emission lines are powered by 
massive stars that collapse out of the cooled gas accumu- 
lated by the cooling flows (Allen 1995). Processes such as 
turbulent mixing layers and shocks in the X-ray absorbing 
medium may also contribute significantly to the observed 
emission- line flux (Crawford & Fabian 1992). 
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Figure 6. A comparison of tlie observed C/D profiles (as sliown in Fig. 2) witii tlie ratios predicted from the deprojection code (Section 
5.2). 



Finally, we note that although the analysis of the B/D 
profiles apparently suggests intrinsic column densities of 
only a few lO'^" atom cm^^ through the centres of cool- 
ing flows, when modelling the absorber as a uniform screen 
in front of the clusters, the true intrinsic column densities 
arc likely to bo 2 — 5 times higher, once systematic effects 
associated due to the spatial distribution of the absorbing 
gas are accounted for (see Section 6). 



5.4 The absorption signature in the 'A' band 

The B/D ratio provides the most reliable probe of the intrin- 
sic column density in the clusters. As noted in Section 1, at 
energies below 0.40 kcV (the lower limit of the B band) the 
point spread function of the PSPC broadens significantly. In 
the Ai (0.1 — 0.2 keV) band a significant fraction of the flux 
is scattered (somewhat asymmetrically) to radii > 1 arcmin. 
In the A2 band, the PSF is better behaved, but the FWHM 
is still ~ 1.7 times broader than in bands B-D. [Note that 
the B/D ratio also has the advantage of being applicable 
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Figure 6. - continued 



to clusters with a wide range of Galactic column densities, 

whereas the A band data arc lost for clusters behind Galac- 
tic column densities > lO'^^ atom cm^^.] 

We have carried out a search for the signature of in- 
trinsic absorption in the A2 data, by forming profiles in the 
ratio A2/T)' , where D' is the D band profile convolved with 
a Gaussian filter of 30 arcsec FWHM. For the two best data 
sets for which this test is applicable, Abell 1795 and Abell 
2199, we determine excess column densities in the central 
bin of ~ 1 — 2 X 10^^ atom cm~^. These column densities 
are significantly less than those inferred from the B/D anal- 



ysis. However such results are to be expected if the absorber 

is spatially distributed throughout the cooling flows, as dis- 
cussed in Section 6, below. Note also that the presence of an 
absorption signature in the A2 data implies that absorbers 
other than Oxygen (which has its K edge at 0.53 keV in low 
temperature gas) must be present in the clusters. 

5.5 The ages of the cooling flows 

Throughout this Section we have limited our comparison 
of the observed and predicted X-ray colour profiles to radii 
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Figure 7. (Upper panels) Comparison of tlie observed C/D colour profiles with the predictions (dotted curves) from the deprojection 
code, as a function of the age of the cooling flows. For Abell 1795, the predictions arc shown for cooling cut-off times of 6.7, 11.4 and 
16.0 X 10^ yr, corresponding to ages of 3.6-6.7 X 10^,6.7-11.4 X 10^, and 11.4-16.0 X 10^ yr, respectively. For Abell 2199 the cooling 
cut-off points arc 5.3, 7.7, 10.3, 12.9 and 15.6 X lO'' yr, implying ages of 3.7 - 5.3, 5.3 - 7.7, 7.7 - 10.3, 10.3 - 12.9 and 12.9 - 15.6 x 10^ 
yr. (Lower panels) The mean cooling time as a function of radius in the clusters. 



where the mean cooling time of the cluster gas is < 5 x 10 
yr. However, for the best data sots, the data from larger 
radii allow constraints to be placed on the age of the cool- 
ing flows (where the age is presumably related to the time 
since the central regions of the clusters were last disrupted 
by a major subcluster merger event). For these data sets, 
the deprojection code has been used to predict C/D profiles 
throughout the clusters, as a function of cooling-flow age. 
At radii where the mean cooling time of the gas is equal 
to or exceeds a specified value, the cluster is assumed to be 
isothermal with no mass deposition occurring. At interior 
radii, where the mean cooling time is less than this value, 
the predictions are calculated in the manner described in 
Section 5.2. 

Fig. 7 shows the observed C/D profiles and the pre- 
dicted results, as a function of cooling flow age, for the two 
best data sets in the sample, Abell 1795 and Abell 2199. 
Also plotted are the mean cooling times as a function of 
radius in the clusters. For Abell 1795, the data arc best de- 
scribed by an age of between 3.6 and 6.7 x lO'"' yr (6.7 x 10^ 
yr being the mean cooling time at the radius where the C/D 
ratio rises, and 3.6 x 10^ yr the mean cooling time for the 
next interior bin.) For Abell 2199, the data suggest an age 
of 5.3 — 7.7 x lO'' yr. The results also imply mass deposition 
rates in the clusters somewhat less than the nominal values 
quoted in Table 3 (which are calculated for an assumed age 
of 1.3 X 10^° yr.) For Abell 1795 we determine a revised mass 



deposition rate of ~ 250 M© yr"^ . For Abell 2199 the re- 
vised mass deposition rate is ~ 80 hlQyr~^ . These results 
are approximately a factor 2 less than the values of 450 
and 160 Mq yr^^ determined under the standard deprojec- 
tion assumptions, but are in good agreement with the rates 
inferred from the analysis of Einstein SSS spectra for these 
clusters by White et al. (1991). 

In general, the data for the other clusters do not place 
firm constraints on the age of the cooling flows. However, for 
Abell 496 the data suggest an age of between 3.3 and 6.8 x 
lO" yr, and a mass deposition rate, M ~ 50 M© yr~^ (again 
in good agreement with the SSS spectral results). For Sersic 
159, the age and mass deposition rate are 6.4 — 10.5 x 10® 
yr and M ~ 220 Mcoyr^^ respectively. For Abell 2597, 
the values arc 2.5 - 5.1 x lO" yr and M ~ 200 M© yr~^ . 
In contrast the data for both Hydra A and Abell 2029 are 
consistent with the cooling flows having been undisturbed 
for close to a Hubble time. 

In detail, the relationship between the age of the cool- 
ing flow and the mean cooling time in the region where 
the C/D ratio changes will be more complex than discussed 
above, and will depend on (undetermined) factors relating 
to the distribution of density inhomogeneities in the cluster 
gas. However, the uncertainties associated with such effects 
should not much exceed those quoted above. 

These results illustrate how standard image deprojec- 
tion analyses, which assume that cooling flows have existed 
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6.1 The effects of partial covering 




Figure 10. The relationship between the apparent excess col- 
umn density, AAf^, and the true intrinsic column density, AA^h, 
determined from the analytical expression in equation 4 (dotted 
curve) and spectral simulations (solid curve). For the spectral 
simulations a Galactic column density of lO^'' atom cm~^ has 
been assumed. Both axes are in units of 10^" atom cm~^. 



for a Hubble time, can overestimate the integrated mass 
deposition rates in clusters by up to a factor of 2. Future 
studies, with the improved spectral and spatial resolution of 
AXAF, will allow such properties to be explored in greater 
detail. The aging of cooling flows may provide a powerful 
probe of the evolutionary history of clusters. 



6 THE INTERPRETATION OF FITTED 
COLUMN DENSITIES 

In this Section we further examine the results on intrinsic ab- 
sorption in the cooling flows. We discuss how measurements 
of intrinsic column densities, determined from fits with uni- 
form foreground screen absorption models, will naturally un- 
derestimate the true intrinsic column densities when the ab- 
sorbing gas is spatially distributed throughout the cooling 
flows. In particular, we explore two limiting distributions 
for the absorbing gas; partial covering (which will apply if 
the absorbing gas is clumped on large scales) and a homo- 
geneous distribution of absorbing material throughout the 
X-ray emitting region, which we refer to as the 'multilayer' 
model. We compare the absorption measurements from the 
ROSAT X-ray colour profile analysis to previous results from 
the Einstein SSS (White et al. 1991) and ASCA. Finally in 
this Section, we comment on the mass of absorbing gas in 
the clusters and on the role of this material as a sink for the 
cooling gas. 



The use of simple uniform-screen absorption models in the 
analysis of cluster spectra can lead to significant underesti- 
mates of intrinsic column densities. A simple illustration of 
this, in the case where the absorbing material covers only a 
fraction, /, of the intrinsic emission, is given in Figs. 8(a) and 
(c). The figures show the column densities inferred from fits 
with simple uniform (foreground) screen absorption models 
to simulated cluster spectra incorporating intrinsic absorp- 
tion with partial covering. For the simulations we assume a 
cluster temperature of 5 keV, a metallicity of 0.5 solar, and 
a redshift, z = 0.01. The cluster is assumed to lie behind a 
Galactic column density, iVn = 10^° atom cm~^, and to be 
intrinsically absorbed by an excess column density, AA^h, 
which covers a fraction, /, of the emitted flux. Intrinsic col- 
umn densities of AA^h = 10^^ and AA^h = 5 x 10^^ atom 
cm~^, and covering fractions of 0.1 — 1.0, were examined. 
[These column densities span the range of values inferred 
from Einstein Observatory SSS and ASCA observations of 
cooling flows; White et al. (1991), Fabian et al. (1994), Allen 
et al. (1996a), Allen et al. (1996b), Fabian et al. (1996).] 

The simple partial covering distribution provides a 
crude approximation to the likely situation in an intrinsi- 
cally absorbed cooling flow. In the simplest case of a geo- 
metrically thin absorbing screen of uniform column density, 
the observed emission can be modelled by two components, 
accounting for the emission from in front of and behind the 
absorber, respectively. If the absorber spans the entire re- 
gion of interest, the components should have equal intrin- 
sic fluxes, and be absorbed by column densities ATh and 
A'h + AA'h. In this case, the covering fraction of the in- 
trinsic absorber is 0.5. If the intrinsic absorber is present 
over only part of the region of interest, as would be the case 
if the extent of the absorber were less than the region being 
analysed, the covering fraction will be less than 0.5. (With 
a real cooling flow, variations in both column density and 
covering fraction are likely.) 

The fits to the simulated spectra were carried out in 
both the 0.1 - 2.0 keV and 0.4 - 2.0 keV energy bands. The 
temperature and metallicity were fixed at 5 keV and 0.5 
Zq, respectively. For covering fractions less than unity, the 
column densities determined from fits with uniform-screen 
(/ = 1) absorption models significantly underestimate the 
'true' absorption. For an intrinsic column density, AA^h = 
5 X 10^^ atom cm~'^, and a covering fraction, / = 0.5, a fit to 
the 0.4 — 2.0 keV spectrum implies a column density of 5.7 x 
10^" atom cm~^ , an excess of only 4.7x 10^° atom cm~^ over 
the Galactic value. Thus the column density inferred from 
the spectral fit to the 0.4—2.0 keV band data underestimates 
the true intrinsic column density by a factor of ~ 10. Even 
more dramatically, when the full 0.1 — 2.0 keV energy range 
is used, the spectral fit implies an excess column density of 
only 4.5 X 10^® atom cm~'^, underestimating the true intrinsic 
column density by a factor of ^ 100 (c/. Section 5.4). In a 
sense, however, the uniform absorption model still provides 
a statistically adequate description of the simulated spectra 
(for signal-to-noise ratios in the simulated spectra typical 
of PSPC cluster observations). Similar results arc obtained 
for the simulations with AA^h = 10^^ atom cm~^, although 
the underestimation of intrinsic column densities from the 
spectral fits is less severe. A fit to the 0.4—2.0 keV spectrum, 
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Figure 8. (a) The apparent column density, N^, as a function of the covering fraction, /, of the intrinsic absorber, determined from fits 
to the simulated spectra with simple uniform-foreground-screen absorption models (Section 6.1). A Galactic column density, Nn = 10^" 
atom cm~^, and an intrinsic column density, ANji = 10^^ atom cm~^, are assumed. We adopt a fixed temperature of 5 keV, a metallicity 
O.SZq and a redshift z = 0.01. Results are shown for fits to both the full 0.1 — 2.0 keV band (lower dashed curve) and the 0.4 — 2.0 keV 

band (upper solid curve), (b) The similar results obtained from the analytical expression in equation 3 for AA^h = 10^^ atom cm^^ and 
photon energies of 0.4 keV (upper solid line) and 0.2 keV (lower dashed line), (c) As for (a) but with AA^h = 5 x 10^^ atom cm~^. (d) 
As for (b) with AATh = 5 x lO^i atom cm-^. 



with / — 0.5, gives a measured column density of 3.9 x 10 
atom cm^"^, implying an excess column density a factor of 
~ 3 less than the true value. 

Such results are not sensitive to the precise emission 
spectrum assumed. If a cooling flow, rather than isothermal, 
model is used both for generating and fitting the simulated 
spectra, results very similar results to those plotted in Fig. 8 
are obtained. Significant discrepancies only axise when dif- 
ferent emission models are used for generating and fitting 
the simulated data. 

In Fig. 9 we show the results from a consistency check 
in which compare the results from a simple, uniform-screen 



fit to a simulated cluster spectrum with A^h ~ 10 atom 
cm^"^, AA^H = 10'^^ atom cm^'^and / = 0.5, with the results 
from a fit to the observed spectrum for the central 2 arcmin 
(radius) region of Abell 1795. (The signal-to-noise ratio in 
the simulated spectrum has been adjusted to approximately 
match that in the Abell 1795 data.) In both cases, the fit to 
the 0.1 — 2.0 keV spectrum implies an excess column density 
of ~ 2 X 10^^ (whereas the B/D results for Abell 1795 imply 
a much larger intrinsic column density of 5 x 10^° atom 
cm-2). 

Some intuitive understanding of these results may be 
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Figure 9. (a) Simulated spectrum of a 5 keV cluster at a redshift z = 0.01, absorbed by a Galactic column density of 10^" atom 
cm~'^ (covering fraction unity) and an intrinsic column density of 10^^ atom cm~^ (covering fraction / = 0.5). The data are shown for 
the full 0.1 — 2.0 keV range, with the best-fitting single temperature, uniform absorbing screen model overlaid. The best-fit column density 
is 1.25j^Q']^^ X lO'^'-' atom cm^'^, implying an excess column density of 2.5^^'^^ X lO-'"'-' atom cm^^. The lower panel shows the residuals 
to the fit. (b) The spectrum for the central 2 arcmin of Abell 1795, at a redshift z = 0.06334, with the best-fitting single temperature, 
uniform absorbing screen model overlaid. The best fit column density is 1.27^g Q4 X 10^" atom cm~^, implying an excess column density 
of l.eto l X lO^S atom cm ^. The measured temperature of 3.04^g'^g keV is lower than the 5 keV assumed in the simulated spectrum, 
due to the strong cooling flow. The Abell 1795 spectrum has the same integrated count rate in the 0.1 — 2.0 keV band as the simulated 
spectrum. Error bars are Icr confidence limits. 



gained from the following argument. We can write the ob- 
served flux as 

where A{E) is the emission at energy E. Then the apparent 
column density determined from a fit assuming a covering 
fraction of unity is 

K = Nj,- log[l - / + /e-''(^'^'^H]/a(i5), (3) 

implying an apparent excess absorption of 

AK = - log[l - / + /e-''(^'^^H]/a(S). (4) 

The absorption cross-section may be crudely approx- 
imated by a{E) « 10-^^(^/0.4) "^/^ cm'^ In Figs. 8(b) 
and (d) we plot N^^ as a function of / for photon energies 
E = 0.20 and 0.40 keV, and intrinsic column densities of 
10^^ and 5 x lO'^^ atom cm^^. (Note that the lowest-energy 
photons included in any fits will dominate the constraints 
on the absorption value.) The analytic curves provide a 
good approximation to the observed behaviour. The appar- 
ent excess column density, AN^, remains low so long as 
significant signal from low energies (below the energy where 
ctNh ~ 1) is included , i.e. E < OAN^^^^i koV, where A^'h = 
10^^ArH,2i cm-2 . Where aANu is small, AN^^ fANn- 

Our results imply that simple fits to PSPC spectra for 
clusters with low Galactic column densities (JVh ~ 10^° 



atom cm"^), which assume covering fractions for the ab- 
sorber of unity, will always return low values for the excess 
column density even if, in reality, / ~ 0.5 and AA^h ^ lO'^^ 
atom cm^'^. If the true intrinsic column densities in clusters 
are > 10"^ atom cm^^, we should only obtain reasonable 
estimates for AA'h with simple spectral fits if we either re- 
strict the bandwidth artificially, or if the cluster has a large 
Galactic column density {Nu ~ 10^^ atom cm"'^). Such a 
scenario occurred naturally in the PSPC studies of Abell 
478 by Allen et al. (1993) and of 2A 0335-^096 by Irwin & 
Sarazin (1995), where significant excess column densities of 
~ 10^^ atom cm"^ were inferred. 

Such arguments also imply a weak positive correlation 
between the apparent excess column density, AN^^, deter- 
mined from fits with uniform foreground screen absorption 
models, and the Galactic column density, A^h. When A^h is 
high then the photon energy dominating the absorption fit 
(generally the lowest energy with a significant detection) is 
also high and cr AA^h tends to be small, so that the inferred 
column density AA^^ ~ fANu- Provided that the covering 
fraction, /, is also high, the model will tend to the 'true' an- 
swer. This is consistent with the results presented in Table 
6, which exhibit a weak correlation between AAT^ and Nn, 
particularly once the radial variation of AA/^ in Abell 426 
is taken into account. 

We have used the relationship between AN^ and AA/h 
in equation 4, and the results from further spectral simu- 
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lations, to estimate the correction factors that should be 
apphed to the ATVh values listed in Table 6, in determin- 
ing the true intrinsic column densities in the clusters. For 
these purposes we assume a covering fraction / = 0.5 for 
the intrinsic absorber (which should approximate the true 
value in the central regions of the clusters). The results are 
shown in Fig. 10. We see that for an apparent excess column 
density, AiV^ = 2 x 10^° atom cm~^, the analytical expres- 
sion and spectral simulations imply true intrinsic column 
densities of 4.5 x 10^" and 6.1 x 10^" atom cm^^, respec- 
tively. Typically, the corrected intrinsic column densities for 
the clusters in our sample are a factor of 2 — 3 larger than 
apparent values determined directly from the B/D analysis. 
[The corrected intrinsic column densities determined from 
equation 4 are listed in Column 4 of Table 6 as AA'h(I).] 

Finally, it is also important to note that the use of 
slightly more complex spectral models incorporating both 
isothermal and cooling-flow components, with the excess ab- 
sorber modelled as a uniform screen in front of the cooling 
flow, will still lead to significant underestimates of the true 
levels of intrinsic absorption in clusters. Typically, only a 
few per cent of the X-ray emission along the line of sight 
to the central 30 arcsec bin will be due to material in front 
of the cooling flows. If the absorbing material is distributed 
throughout the cooling flow, as required by the observations 
reported here, any model which places the absorber entirely 
in front of the cooling flow will underestimate the true in- 
trinsic column densities in much the same manner as in the 
simulations discussed above. 

6.2 The multilayer absorption model 

Wc have also explored a second plausible distribution for 
the X-ray absorbing material in which the absorbing gas is 
assumed to reside in many small clouds, each of which has an 
individual column density much less than the total intrinsic 
column density, AA^h. The clouds are assumed to be evenly 
distributed throughout the X-ray emitting plasma. For an 
intrinsic emission spectrum A{E), the observed spectrum, 
A'{E), emerging from the emitting/absorbing region may 
be written as 

A series of spectra incorporating absorption as de- 
scribed by equation 5 (which we hereafter refer to as the 
'multilayer' absorption model) were simulated. Again, the 
intrinsic emission spectrum was taken to be a 5 keV ther- 
mal plasma, with a mctallicity, Z — 0.5Zq, and a redshift, 
z = 0.01. Total intrinsic column densities of AA^h = 10^^ 
and 5 x 10^^ atom cm~^were examined. 

The simulated spectra were then fitted with simple uni- 
form foreground screen absorption models. As with the par- 
tial covering model discussed in Section 6.1, the apparent 
excess column densities determined with the uniform ab- 
sorbing screen model significantly underestimate the true 
intrinsic column density. For a Galactic column density, 
ATh = 10^° atom cm , and intrinsic multilayer column den- 
sities of AA^H = 10^"*^ and 5 x 10^^ atom cm^^, the apparent 
excess column densities determined with the uniform ab- 
sorbing screen model are AAT^ = 2.65 and 10.3 x 10^° atom 
cm~^, respectively. Thus, with a multilayer distribution for 



the intrinsic absorbing gas, the simple uniform foreground 
screen model underestimates the true intrinsic column den- 
sity by factors of ~ 4 and 5 respectively. Note also that 
these factors are essentially independent of the Galactic col- 
umn density, since multilayer absorption, with total column 
densities of AA'h = 10^^ and 5 x 10^^ atom cm~^, can be 
well-approximated by uniform absorbing screens with col- 
umn densities of ~ 0.26 and 1.0 x lO^'^ atom cm~^, respec- 
tively. In Column 5 of Table 6 [AA^h(2)], we list the intrinsic 
column densities for the sample of cooling flows in this pa- 
per, corrected (with an intermediate correction factor of 4.5) 
for an assumed multilayer distribution of the absorbing gas. 

6.3 The comparison of results from ROSAT, SSS 
and ASCA 

The results presented here confirm that intrinsic X-ray ab- 
sorption is a common property of cooling flows. All of the 
clusters with previous positive detections of intrinsic absorp- 
tion, from the Einstein SSS study of White et al. (1991), also 
exhibit absorption signatures in their ROSAT X-ray colour 
profiles. The column densities inferred from the two studies 
are also similar, ranging from a few 10^° to a few lO'^^ atom 
cm~^. However, although the results reported here (Table 
6 columns 4, 5) include corrections for the effects of the 
likely spatial distribution of absorbing gas on the measured 
column densities, the White et al. (1991) results were all 
determined under the more simple assumption that the ab- 
sorber acts as a uniform screen in front of the clusters, with 
a covering fraction of unity. When the same simple model 
for the distribution of the absorbing gas is used with the 
ROSAT data, (Table 6, Column 3) we obtain results that 
are typically factors of 2 — 6 less than the SSS results. In 
this Section we examine whether this discrepancy should 
simply be attributed to calibration uncertainties, or if it in- 
stead provides further information about the nature of the 
absorbing gas. 

The SSS observations were made in the 0.6 — 4.5 keV 
energy band, with a spectral resolution of ~ 160 eV, a higher 
energy band and a higher spectral resolution than that of 
the PSPC. We have examined the variation in the results 
from fits to simulated cluster spectra, incorporating intrinsic 
absorption with both partial covering and multilayer distri- 
butions, when observed with the two instruments. (With the 
SSS simulations, the effects of ice build-up on the instrument 
response were accounted for.) Total multilayer and partial 
covering column densities of lO^'^ — 10^^ atom cm~^ were 
examined, with a default Galactic column density of 10^° 
atom cm~^ assumed. The simulated spectra were then fit- 
ted with simple uniform foreground screen models, assuming 
a covering fraction for the excess absorption of unity. 

With a multilayer distribution for the intrinsic absorb- 
ing gas in the simulated spectra, the fits with the uniform- 
screen model give very similar results for the excess column 
densities with the PSPC and SSS instruments (in both cases 
the fit results underestimate the true intrinsic column den- 
sities by factors of ~ 4 — 5; cf. Section 6.2). However, when 
the intrinsic absorber in the simulated spectra is distributed 
according to the partial covering model, clear differences be- 
tween the PSPC and SSS results become apparent. 

The differences between the PSPC and SSS results, 
from the fits to the simulated spectra, incorporating partial 
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Figure 11. The ratio of the apparent excess column densities, 

AA^^, inferred from fits witli simple uniform foreground screen 
absorption models to simulated SSS and PSPC spectra incor- 
porating intrinsic absorption with a covering fraction, /, of be- 
tween 0.1 and 1.0. The solid curves show the results (from top 
to bottom) for AA^h = 10,5 and 1 X 10'^-'^ atom cm~^, with a 
fixed Galactic column density, TVh, of 10-^° atom cm~^. The dot- 
dashed curve shows the results for AATh = 10^^ atom cm~^ and 
Nil = 3 X 10^" atom cm~^. The dotted curve shows the ratio of 
the results for the simulated ASCA SIS and PSPC spectra, for 
AA^H = 5 X 10^1 atom cm"^ and A^h = 10^'' atom cm^^ ^o^g 
the agreement between the SIS/PSPC and SSS/PSPC results. 

covering, are summarized in Fig. 11. For a Galactic column 
density of ~ 10^° atom cm~^, and intrinsic absorption with 
a column density of a few lO^'^ atom cm~^ and a covering 

fraction / = 0.2 — 0.6, discrepancies of factors of between 2 
and 6 will naturally arise between measurements made with 
the SSS and PSPC instruments (in the sense the SSS results 
exceed the PSPC results). Such an effect is sufficient to ac- 
count for the bulk of the discrepancy between the Wfiitc et 
al. (1991) SSS results and the ROSAT results reported here. 

We have also carried out a further comparison, between 
the PSPC results and the results from simulated spectra 
generated with an instrument response appropriate for the 
Solid-state Imaging Spectrometers (SISs) on ASCA. The 
dotted curve in Fig. 11 shows the ratio of the excess col- 
umn density determined from the ASCA SIS simulations 
to that with ROSAT PSPC, for a Galactic column density 
iVn = 10^° atom cm and an intrinsic column density of 
5 X 10^^ atom cm~^. The results are very similar to those 
obtained from the comparison of the SSS and PSPC data, 
for the same absorption parameters. 

Since the partial covering model can account for many 
of the observed properties of the intrinsic absorber, we have 
made one further check on the validity of this model by 
applying it to observed ROSAT PSPC, Einstein SSS and 
ASCA SIS spectra for Abell 2199 (which has the best com- 



bined data set of any of the clusters in the present sample). 
Spectra from the central 2 arcmin (radius) region of Abell 
2199 were extracted from the PSPC and ASCA SIS data 
sets (2 arcmin radius being the smallest region that can be 
consistently analysed with the two instruments, given the 
relevant PSFs). The SSS offers no spatial resolution but did 
have a 3 arcmin (radius) circular aperture, reasonably well- 
matched to the regions studied with the other instruments. 
With the ASCA analysis, the data from both SIS detectors 
were used. With the SSS analysis, the data from all three 
observations of the source (White et al. 1991) were included. 

The spectra were fitted with a model incorporating 
emission from an isothermal plasma and a cooling flow of 
80 M0yr"^ (c/. Section 5.5), at a redshift oi z = 0.0309, 
absorbed by a (fixed) Galactic column density of 8.7 x 10^^ 
atom cm~^ (except for the PSPC data where the Galactic 
column density was fixed at 5 x 10^^ atom cm~^; see Section 
5.3) . A second absorption component of variable column 
density and covering fraction was also assumed to act on 
the emission. The cooling flow was forced to cool from the 
ambient cluster temperature, which was a free parameter in 
the flts (with the exception of the ROSAT analysis where, 
due to the limitations of the data, the upper temperature of 
the flow was flxed at 3.5 keV, in agreement with the results 
of Section 5). The metallicity of the X-ray emitting gas and 
the emission measure of the isothermal component were also 
free parameters in the flts. With the ROSAT data, the spec- 
tral analysis was limited to the 0.4 — 2.0 keV energy range. 
For the ASCA SIS and Einstein SSS data, spectral ranges of 
0.6 — 10.0 keV and 0.6 — 4.0 keV were used, respectively. All 
appropriate ice corrections were applied to the SSS data. 

The joint confidence contours on the column density 
and covering fraction of tlK> excess absorber in Abell 2199 
are shown in Fig. 12. The contours indicate the regions of 
68, 90 and 99 per cent confidence on two interesting pa- 
rameters. The SIS and SSS data, in particular, exhibit good 
agreement. For the SIS data, the best-fit excess column den- 
sity is 2.5lj^ J X lO'^^ atom cm^'^ and the covering fraction, 
/ = 0.53lo For the PSPC data, the best-fit values are 
AATh = 1.2+^:^ X 10^^ atom cm"^ and / = 0.431^:^3, where 
errors are the 90 per cent (Ax'^ = 2.71) confidence limits on 
a single interesting parameter. 

The simple partial covering model thus provides a rear- 
sonable description of the excess absorption in Abell 2199 
and leads to reasonable agreement in the parameter values 
inferred from the different data sets. However, although the 
errors on the individual parameter values overlap, the Ax^ 
contours shown in Fig. 12 indicate an inconsistency between 
the PSPC data and the other two data sets. This implies ei- 
ther a deficiency in the spectral model or a calibration prob- 
lem with the instruments. The simple partial-covering model 
used undoubtedly underestimates the complexity of the ab- 
sorption in a real cluster, where the column density acting on 
the emission from different regions is likely to span a range of 
values. The abundances of separate elements in the absorber 
could also deviate from the assumed solar values, perhaps 
as a function of radius. Further complexities would be in- 
troduced if the absorbing material were dusty, which recent 
evidence suggests is likely to be the case (Fabian, .Johnstone 
& Daines, Allen et al. 1995, Voit & Donahue 1995, Allen et 
al. 1996). If, however, the origin of the discrepancy is a cal- 
ibration problem, the agreement between the SIS and SSS 
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Figure 12. The joint confidence contours on the column density and covering fraction of the intrinsic absorber in Abell 2199 for the 
ASCA SIS (soUd curves), ROSAT PSPC (dashed curves) and Einstein SSS (dotted curves) data. The contours mark the Ax^ = 2.30, 4.61 
and 9.21 limits, corresponding to 68, 90 and 99 per cent confidence for two interesting parameters. 



results would suggest that the problem is likely to lie with 
the PSPC data. 



6.4 The mass of absorbing gas 

The results on the intrinsic column densities, summarized in 
Table 6, allow a simple estimation of the mass of absorbing 
gas in the central 30 arcsec (radius) regions of the clusters 
to be made. Adopting the partial covering model for the 
distribution of the absorbing gas, with / = 0.5, and assum- 
ing solar metallicity in this material, the mass of gas can be 
written as M^bs = 1-62 x 10*'r^AiVH,2o M©, where r is the 
radius of the 30 arcsec region in kpc and AN-h,20 is the in- 
trinsic column density in units of lO'^" atom cm~^. We can 
also estimate the time, tacc, required for the cooling flows 
to accumulate these masses of gas, given the observed mass 
deposition rates in these regions (tacc = Aiabs/M). [Note 
that the mass of absorbing gas is actually a projected mass 
through the central 30 arcsec region, whereas the mass de- 
position rates apply to the central 30 arcsec volume. Thus, 
where significant mass deposition occurs outside the central 
30 arcsec region, the accumulation timescales, tacc, quoted 
will slightly overestimate the true values. Assuming M oc r, 
the projected mass deposition rate through the central 30 
arcsec region will bo ~ 30 per cent higher than the rate 
within that volume, and the accumulation tirncscale will be 
~ 25 per cent shorter than the quoted result.] The masses 
of absorbing gas, mass deposition rates, accumulation times, 
and the masses of hot X-ray emitting gas within the central 
30 arcsec (radius) volumes, are summarized in Table 7. 

Prom Table 7 we see that the accumulation timescales 
are typically only a few 10* years, much less than the ages of 



the cooling flows (c/. Section 5.5.) This implies that most of 
the mass deposited by the cooling flows in the central regions 

of the clusters must reside in some form other than X-ray ab- 
sorbing gas (perhaps low mass stars or brown dwarfs; Fabian 
1994). Note, however, that if the metallicity of the absorb- 
ing material is substantially sub-solar, the true masses of 
absorbing gas are likely to be ~ 2 — 3 times larger than the 
values listed in Table 7. Note also that the mass of absorb- 
ing material (assuming solar metallicity) is typically < 10 
per cent of the mass of the X-ray emitting gas in the same 
regions (although for Abell 478 and MKW3s the values are 
somewhat higher.) 



7 CONCLUSIONS 

In this paper we have described and discussed the results 
from a detailed X-ray colour proflle analysis of a sam- 
ple of nearby, luminous cooling flow clusters observed with 
the ROSAT PSPC. Our analysis of the C/D (0.80 - 1.39 
keV/1.40 — 2.00 keV) ratios clearly demonstrates the pres- 
ence of distributed cool gas in the central regions (r < rcooi) 
of the cooling flow clusters. The spatial distributions and 
emissivity of the cooling gas arc in excellent agreement with 
the predictions from standard image deprojection analyses. 
In the outer regions of the clusters the data suggest approxi- 
mate isothermality. Our results place new constraints on the 
ages of the cooling flows and provide the most accurate de- 
terminations of the mass deposition rates in these systems 
of any imaging study to date. 

The analysis of the B/D (0.41 - 0.79 keV/1.40 - 2.00 
keV) ratios provides flrm evidence for spatially distributed 
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Figure 13. A summary of the results on the distributions of eool gas and intrinsie absorbing material in Abell 426, Abell 1795 and 
Abell 2199. The y axes show the C/D and B/D ratios normalized to the values at large radii (Table 5). The x axes are in units of the 
cooling radii, r^^^i (see Table 3). The C/D profiles clearly demonstrate the presence of distributed cool gas within the cooling radii, and 
the B/D data the need for excess absorption in the central regions. Other details as in Figs. 3 and 5. 
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Table 7. The masses of absorbing gas and accumulation timescales 



Cluster 




M 


^acc 


Mhot 




(109 Mo) 


(Mgyr-i) 


(10* yr) 


(iQio Mo) 


Abell 85 


15.0+tI 

— 4.4 




4.1+1-3 


20.9 ± 0.4 


Abell 426 


— 0.4 


16.0+?:9 


0.25+0-^1 

— 0.25 


3.0 ±0.1 


Abell 478 


143+19 


160+^ 


8.9+{i 


88.6 ±1.0 


Abell 496 




1' -^-0.7 


4.0+{-l 


6.8 ±0.1 


Hydra A 


21.11-9 


69.4+1-^ 


n Q+0.6 


25.9 ±0.3 


Abell 1795 


20.012;^ 


83.4+2-0 


2 4+0-3 

^•^-0.2 


35.7 ±0.3 


Abell 2029 






0.88l0-|^ 


57.7 ±0.9 


MKW3s 


34.7+^-1 


I9.7t\t 


17 6+2-9 


9.0 ± 0.3 


Abell 2199 


"•'-0.7 


13.6l«:^ 


4.9+0-? 

— O.D 


5.6 ± 0.1 


Cygnus A 


44+11 


lists 


q 7+0.9 
•^•'-0.8 


39.3 ±0.6 


Sersic 159 


1 Q Q+3.8 
10.0_3 g 


61.6l^-g 


9 9+0.6 


24.6 ±0.3 


Abell 2597 


oc:+10 
•3O_10 


126+^ 


2 8+0-* 


72.5 ±1.1 



The masses of absorbing gas, Af^bs (calculated assuming solar mctallicity in the absorbing material), mass deposition rates, M, accumu- 
lation timescales, tacc, and masses of hot, X-ray emitting gas, Afbot; in the central 30 arcsec regions of the clusters. The Af^bs results are 
derived from the [AAfH(l)] column density measurements made with the partial covering model (Table 6, Column 4). A covering fraction 
/ = 0.5 has been assumed for all clusters except MKW3s, for -which / = 0.75 was used. The M and Mhot values apply to the central 30 
arcsec (radius) volume. 



intrinsic absorbing material in the cooling flows. The absorb- 
ing gas generally exhibits an increasing column density with 
decreasing radius, and is confined to radii < r-cooi- Our re- 
sults arc consistent with the largo intrinsic column densities 
(AA^H ~ 10^1 atom cm ) inferred from previous studies 
with the Einstein Observatory SSS and ASCA (White et al. 
1991; Fabian et al. 1994) but require that the absorber only 
partially covers the emitting region. This suggests significant 
clumping of the absorber on large scales. 

We have discussed the complexities involved in the mea- 
surement of intrinsic column densities in clusters and have 
shown that the use of simple spectral models, which treat 
the absorber as a uniform foreground screen (in front of the 
cooling flow or cluster) will naturally lead to significant un- 
derestimates (by factors of ~ 2 — 5) of the true amounts of 
absorbing material. 

Assuming solar metallicity in the absorbing material, 

the masses of absorbing gas in the central 30 arcsec regions of 
the clusters typically correspond to < 10 per cent of the mass 
of the hot. X-ray emitting gas. The timescales required for 
the cooling flows to accumulate such masses are significantly 
shorter than (~ 10 per cent of) the likely ages of the cooling 
flows, and imply that the bulk of the matter deposited by the 
cooling flows (in the central 30 arcsec radius regions) resides 
in some form other than X-ray absorbing gas, presumably 
stars or sub-stellar objects. 

The new results presented in this paper provide strong 
support for the standard inhomogeneous model of cooling 
flows in clusters of galaxies. Future studies of the range of 
temperature/density phases in cooling flows (as a function 
of radius), more precise determinations of the 'ages' of cool- 
ing flows, and more detailed mapping of the intrinsic col- 
umn densities in clusters, will become possible following the 
launch of AXAF. 
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